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A new semiempirical formula, with only three parameters, 
is proposed for cluster decay half-lives. The parameters of 
the formula are obtained by making a least squares fit to 
the available experimental data. The calculated half-lives are 
compared with an earlier proposed model-independent scaling 
law. Also, the calculated results of this formula are compared 
with the recent results of the preformed cluster model for 12 C 
and 14 C emissions from different deformed and superdeformed 
Nd and Gd parents. The results are in good agreement with 
experiments as well as other models. 
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I. INTRODUCTION 

In a radioactive decay series, the end product is 
reached not only via the emission of a and (3 particles 
but also directly via heavy nuclei emission, the clusters, 
like Carbon, Oxygen, Florine, Neon, Magnesium and Sil- 
icon. Such a process is known as cluster decay, first pro- 
posed theoretically in 1980 by Sandulcscu, Poenaru and 
Greiner [1], before its experimental realization in 1984 
by Rose and Jones [2]. In this decay process a parent 
nucleus (A,Z), with mass number A and charge num- 
ber Z, breaks into two fragments, viz., the emitted clus- 
ter (A 2 ,Z 2 ) and the associated daughter (Ai,Zi), where 
A=Ai+A 2 , Z=Zi+Z 2 . The light fragment (A 2 ,Z 2 ) is a 
cluster, heavier than the a-particle but lighter than the 
lightest fission fragment observed so far. 

The detection of cluster-decay was hindered mainly 
due to large pile-up of a particles, but with improved 
facilities these difficulties are removed and a large num- 
ber of cluster-decays are observed from different radioac- 
tive nuclei. In the years that followed, the cluster decay 
process has been studied extensively using different theo- 
retical models with different realistic nuclear interaction 
potentials. In general, two kinds of models are used for 
explaining the observed and / or for predicting new decay 
modes. In one kind of these models, the a-particle as well 
as the heavy cluster (s) were assumed to be pre-born in 
a parent nucleus, before they could penetrate the barrier 
with the available Q-valuc. These models are called the 
preformed cluster models [3-7] . In such a model, the clus- 
ters of different sizes (mass and/ or charge numbers) are 
considered to be preformed in the parent nucleus, with 



different probabilities. The Gamow-like barrier penetra- 
tion is also taken into account in these models. In the 
other kind of models, only the Gamow's idea of barrier 
penetration is used, without considering the cluster(s) 
being or not being preformed in the parent nucleus. In 
other words, in this kind of models, called the unified 
fission models [1,8-13], the cluster radioactivity is con- 
sidered as a simple barrier penetration phenomenon, in 
between the a-decay and the spontaneous fission. In this 
paper, we attempt to give a model- independent, semiem- 
pirical formula for studying this above mentioned process 
of exotic cluster decay. 

Geiger and Nuttal [14] were the first who proposed a 
semiempirical law connecting the a-decay half-life and 
its Q-value. Now, a large number of a-decays are ob- 
served from medium mass to superheavy nuclei, and sev- 
eral attempts have been made to give a universal formula 
[15-17]. These formulae vary among themselves mainly 
in the number of parameters. One such scaling law, pro- 
posed recently by Horoi et al [15], accounts for both the 
a and cluster decays, and is given as, 



logT 1/2 = (ai/i x + &i) 



- 7 



+ {a 2f i x + b 2 ). 



(1) 



This scaling law involves six parameters a± = 9.1, b\ = 
-10.2, a 2 = 7.39, b 2 = -23.2, x = 0.416 and y = 0.613, 
which were obtained by fitting 119 a-decays and 11 
cluster-decays from various even-even parents. In the 
following, we propose a new semiempirical formula based 
on only three parameters that are least squares fitted to 
the cluster-decay data alone. So far, no attempt has been 
made to fit the a-decay data. 

The evolution of our proposed formula is presented 
in Section II, followed by a brief description of the pre- 
formed cluster model (PCM) of Gupta and collaborators 
[3-5] in Section III. The PCM is recently used by two of 
us and collaborators [18] to calculate the cluster-decay 
half-lives of some deformed and superdeformed Gd and 
Nd parents. Some of these results of Ref. [18] are used 
here for comparisons with the calculations based on the 
proposed new semiempirical formula. The results of our 
calculation are discussed in Section IV, and the summary 
and conclusions are presented in Section V. 



1 



II. THE NEW SEMIEMPIRICAL FORMULA 



III. THE PREFORMED CLUSTER MODEL 



We base the new formula for cluster decay half-lives on 
the following three simple experimental facts: 

(i) It is known from experiments that the cluster de- 

cay half-lives increase with the size (mass and/or 
charge) of the clusters. Hence, the empirical for- 
mula should contain terms showing direct depen- 
dence on the mass number and charge number of 
the cluster. 

(ii) The same cluster is emitted by different parents and 

hence the formula should contain dependence on 
the mass and charge asymmetries 



7] = 



respectively. 



A x - A 2 



Vz = 



Z\ — Z 2 



(2) 



(iii) Since the a-decay and cluster-decay are physically 
similar processes, the Q-dependence is taken to be 
the same as in Geiger-Nuttal law for a-decay, i.e., 
logT 1/2 ex Q-V2. 

Combining the above three results, we get 



logT^ = 



AZ _ aA 2V + bZ 2 rj z 



+ c, 



(3) 



where the constants a — 10.603, b = 78.027 and c = 
—80.669 are obtained by making a least squares fit of 
the available experimental half-lives for exotic cluster 
decays alone, with an rms deviation rf rTOS =0.89(s), de- 



i=l 



with / (x i , dj ) de- 



fined as d rm 

noting the function in Eq. (3), n the number of measure- 
ments and yi the experimentally observed values. The 
a i = ~^27=i(yi ~ y) 2 ' the var i anc e, gives the standard 
deviation a%, with y = — Y^j—i Vi giving the arithmetic 
average of the expermentally measured quantities. We 
refer to this formula as the AZ-formula {AZF). Note that 
the dependence on both r] and r\z must be included in the 
formula (3) since these are separately measurable quan- 
tities, though the coupling between them is known to be 
weak [19-21], as is also evident when we consider the 77 
and r]z dependence separately. These special cases of AZ- 
formula, i.e., 6=0 or a=0 in (3), respectively, are expected 
to give reasonably good results, though poorer than the 
AZF results. These truncated expressions are referred to 
as the AF and ZF formulae in the following, whose con- 
stants arc also obtained directly by the least squares fit 
to data. These are: a = 30.568 and c = -51.348 for AF 
and b = 112.197 and c = -89.025 for ZF, with rms devia- 
tions d rms =l. 652 and 1.112 (s), respectively. Apparently, 
the rms deviations for truncated expressions are larger, 
and hence the fits are poorer, compared to the total ex- 
pression (3). This will also be evident when we compare 
the results of these expressions with experimental data 
in Section IV. 



In the prefomed cluster model (PCM) of Gupta and 
collaborators [3-5], the decay half- life is given in terms 
of three factors, the Po, P and v, as 



logT 1/2 = ln2/P Pu, 



(4) 



where Po is the preformation probability of the two frag- 
ments (the cluster and daughter nuclei) in their respec- 
tive ground states, P the probability to tunnel the con- 
fining nuclear interaction barrier and v as an assault fre- 
quency. 

For calculating Po and P, Gupta et al. introduced 
the coupled motion in dynamical collective coordinates 
of mass asymmetry r\ and relative separation R via a sta- 
tionary Schrodinger equation 



H{n,R)^{ri,R)=Eil>{ri,R), 



(5) 



with the potential in it defined by the sum of experimen- 
tal binding energies [22], the Coulomb and the proximity 
[23] potentials, as 

V{7 1 ,R) = Y,B i {A i ,Z i ) + ^^+V P . (6) 

i—l 

Here the charges Z^ are fixed by minimizing the potential 
(without Vp) in the charge asymmetry co-ordinate ij z . 
Equation (5) is solved in the decoupled approximation, 
which gives 



Po oc |V(r?)| 5 



Pcx |V(P)| 2 



(7) 



(8) 



Only the ground state solution is relevant for the clus- 
ter decay to occur in the ground state of the daughter 
nucleus. Then, for ^-motion, the properly normalized 
fractional preformation probability for, say, cluster A 2 
at a fixed R (=R a =C t =Ci+C2, Ci being the Sussmann 
central radii) is 



P (A 2 ) = \i,( v )f 



(9) 



For R-motion, the potential V(R) is obtained from Eq. 
(6) for fixed rj and, instead of solving the corresponding 
radial Schrodinger equation, the WKB approximation is 
used for calculating the penetrability P. Finally, the as- 
sault frequency v in PCM is defined by considering that 
the total kinetic energy, shared between the two frag- 
ments, is the positive Q-value: 



v = v /Pq 



y/2Q/mA 2 
Po 



(10) 



Here Po is the spherical radius of the parent nucleus and 
m, the nucleon mass. 
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IV. RESULTS AND DISCUSSION 

Figure 1 and Table I give our calculated logarithms of 
decay half-lives by using the AZ-formula (3) for different 
clusters emitted from various radioactive parents, com- 
pared with the experimental data. In Fig. 1, we have 
also plotted the results of our calucations for AF (b = 0) 
and ZF (a = 0) versions of the AZ-formula. It is evident 
that the AZF fit to the data is better than the AF and 
ZF fits, as expected from our discussion above in Section 
II. 

In Table I, we have also added the experimental data 
on Q- values and the results of another calculation using 
Eq. (1) of Horio et al. [15]. The comparison between the 
experiments and formulae (1) and (3) are also displayed 
in Fig. 2 for the illustrative cases of 14 C and 24 Ne cluster 
decays (respectively, the upper and lower panel) . Appar- 
ently, our semiempirical AZ-formula is much closer to 
experiments, as compared to the other formula due to 
Horio et al. [15]. 

Finally, in order to compare the results of our semiem- 
pirical formula with the results of a model-dependent the- 
oretical cluster decay calculation, we use a recent calcu- 
lation [18] of two of us (MB and RKG) and collaborators 
based on PCM for the emission of 12 C and 14 C clusters 
from ^ 3 ~ 137 Nd and ^~ 158 Gd parents. It may be men- 
tioned here that these cluster-decay calculations are of 
interest only for nuclear structure information since their 
experimental observation may not be feasible in the near 
future [18]. Figs. 3(a) to (d) show the results of the PCM 
calculations, compared with those from the semiempiri- 
cal AZ-formula (3) and the scaling law (1) of Horoi et 
al. [15]. We notice that, in general, the predictions of 
our semiempirical law lie higher than those of the PCM 
and Horoi et al. This is more so for Nd parents than 
for Gd parents where the predictions of the three calcu- 
lations are nearly similar. The interesting point is that 
the three calculations predict an exactly the same struc- 
ture for logTx/2 vs. A, the parent mass number. For 
quantitative comparisons, it may be reminded here that 
the constants a, b and c of our scmicmprical formula are 
obtained by fitting the data from trans-actinide region, 
which may or may not be good for this trans-tin region. 



good comparisons suggest that this formula could be used 
to make predictions for the guide of new experiments, 
similar to that of Poenaru et al. [41]. 
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Figure Captions 

Fig. 1 The logarithms of decay half-lives for different 
clusters emitted from various radioactive parents, 
calculated by using AZ-formula (AZF) and com- 
pared with experimental data. Also, the results 
of calculations for AF (b=0) and ZF (a=0) trun- 
cations of AZF are shown for comparisons. The 
parents are labelled at the top X-axis and the cor- 
responding clusters emitted by these parents are 
labelled at the bottom X-axis. 

Fig. 2 The experimental data on logarithms of decay 
half-lives for the emission of 14 C and 24 Ne clusters 
from different radioactive parents, compared with 
the results of calculations using AZ-formula (AZF) 
proposed here and the scaling law of Ref. [15]. 

Fig. 3 The logarithms of half-lives for the emission of 
12,1 C clusters from different deformed and su- 
perdeformed Nd and Gd parents. Our calculations, 
using AZ-formula (AZF), are compared with the 
calculations based on PCM and the scaling law of 
Horoi et al. [15]. 



V. SUMMARY 



In this paper, we have proposed a model-independent 
three parameter formula for calculating the half-lives of 
cluster decays of nuclei. The evolution of the formula is 
based on three simple experimental observations about 
the characteristics of exotic cluster decays. The inputs 
of the formula are simply the mass and charge numbers of 
the parent and cluster, along with the Q-value of decay. 
The predictions of the proposed formula are comparable 
to another model-independent scaling law as well as to 
the well known preformed cluster model. The resulting 
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TABLE I. The logarithms of decay half-lives of different 
clusters emitted from various radioactive nuclei, calculated by 
using the semiempirical AZ-formula (3) and the scaling law 
(f) compared with the experimental data. The experimental 
Q-values are also given and the data are tabulated in the 
increasing order of the mass number of the clusters. The last 
column gives the reference to the source of experimental data. 



Parent Cluster Daughter Q 


Expt - (MeV) log 


'.T$ F (s) 


log T*f [16 l(s) 


log T*T' (s) 


Ref. 


221 Fr 




± I 


31.28 


14.54 


13.56 


14.52 


[24] 


221 Ra 




207 Pb 


32.39 


12.96 


12.28 


13.39 


[24] 


222 Ra 




208 Pb 


33.05 


11.98 


11.00 


11.01 


[25] 


223 Ra 


14 


209 Pb 


31.85 


13.81 


13.38 


15.20 


[25] 


224 Ra 

Till 




210 p. 


30.54 


15.94 


16.13 


15.68 


[26] 


225 Ac 


14^» 


211 Si 


30.48 


16.20 


17.26 


17.16 


[27] 


226 Ra 


14 c 


212 P6 


28.21 


20.08 


21.50 


21.19 


[28] 


22Srpr 
lit 


20 n 


208 p ft 


44.72 


21.90 


21.20 


20.72 


[29] 


230 jj 


22 Ne 


208p fe 


61.59 


21.78 


19.28 


20.14 


[30] 


231 Pa 
I a 


23 p 


208 p, 

I U 


51.84 


24.30 


23.85 


26.02 


[31] 


230 Th 


24 Ne 


206 Ha 


57.78 


25.77 


23.88 


24.61 


[32] 


231 Pa 
i a 


24 Ne 


207 Tl 


60.42 


23.62 


21.30 


23.23 


[33] 


232 jj 


24 iVe 


208 Pb 


62.31 


22.24 


19.94 


21.08 


[34] 


233 jj 


24 iVe 


209 Pb 


60.5 


23.87 


22.53 


24.83 


[32] 


234 u 


24 iVe 


2w Pb 


58.84 


25.42 


25.01 


25.92 


[35,36] 


235 jj 


24 iVe 


211 Pb 


57.36 


26.87 


27.31 


27.42 


[35,37] 


233 jj 


25 Ne 


20Sp b 


60.75 


24.15 


22.98 


24.83 


[32] 


235 jj 


25 Ne 


210 Pb 


57.83 


26.94 


27.49 


27.42 


[35,37] 


234 u 


26 Ne 


208 Pb 


59.47 


25.85 


25.75 


25.92 


[35,36] 


234 u 


2S Mg 


206 Hg 


74.13 


26.24 


24.74 


27.54 


[37] 


236 Pu 


2S Mg 


20S Pb 


79.67 


23.01 


20.83 


21.67 


[38] 


236 jj 


2S Mg 


208 Hg 


71.69 


28.18 


27.98 


27.58 


[39] 


238 Pu 


2S Mg 


2W Pb 


75.93 


25.70 


25.39 


25.70 


[40] 


236 jj 


30 Mg 


20e Hg 


72.51 


28.36 


28.36 


27.58 


[39] 


238 Pu 


30 Mg 


208 Pb 


77.03 


25.71 


25.41 


25.70 


[40] 


238 Pu 


32 Si 


206 Hg 


91.21 


25.99 


25.68 


25.27 


[40] 
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Parents 



CM 



CD 
O 



35 



22V 222 -r» 224 -r» 22 6-r» 230, t 230™ 232 T T 234 T T 233 T T 234 T T 236,, 238,, 238,, 

Fr Ra Ra Ra U Tn U U U U Pu Pu Pu 
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25 - 



20 - 



15 - 
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20.- 23-,-, 24, T 24, T 24, T 25, T 28,, 28,, 30,, 32 n . 

j F | Ne, Ne, Ne, Ne, Mg | Mg | Mg | S i 

14 14 14 14 22, T 24, T 24, T 24, T 25, T 26, T 28, , 28, , 30, , 

C C C C Ne Ne Ne Ne Ne Ne Mg Mg Mg 

Clusters 

Fig.l "A new semiempirical formula for..." by M. Balasubramaniam et al. 



Fig.2 "A new semiempirical formula for..." by M. Balasubramaniam et al. 




Fig. 3 "A new semiempirical formula for..." by M. Balasubramaniam et al. 




